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Abstract 

The present theoretical status of neutrino interactions in the few-GeV region is reviewed. Quasielastic scattering, 
pion production, photon emission and their importance for neutrino oscillation studies are discussed, making emphasis 
on the open questions that arise in the comparison with new experimental data. 

Keywords: neutrino-nucleus reactions, form factors, baryon resonances, quasielastic scattering, pion production, 
photon emission 



1. Introduction and motivation 

Recent years have witnessed an intense experimental 
and theoretical activity aimed at a better understanding 
of neutrino interactions with nucleons and nuclei. While 
the main motivation for these efforts is the demand from 
oscillation experiments in their quest for a precise de- 
termination of neutrino properties, the relevance of neu- 
trino interactions with matter is more far-reaching. They 
are important for astrophysics, physics beyond the stan- 
dard model, hadronic and nuclear physics. 

In the few-GeV neutrino-energy region, where most 
oscillation experiments operate, the dominant reaction 
channel through which neutrinos reveal themselves (and 
their flavor) is charged-current quasielastic scattering 
(CCQE) V; n — > r p. Oscillation probabilities depend 
on the neutrino energy, unknown for broad fluxes and 
usually obtained from the measured angle and energy 
of the outgoing lepton using two-body kinematics 



2 [nin -Ei+ ^Ef - cos 6»,j 

This determination is only exact for free neutrons and 
under the condition that inelastic events (mainly pion 
production ones) are identified and excluded. For nu- 
clear targets, as in all modern neutrino experiments, the 
situation is more involved as shown in Fig.[T| Here the 



reconstructed neutrino-energy distribution is obtained 
applying Eq. ([T]i to charged-current (CC) events gener- 
ated with the GENIE Monte Carlo (MC) generator U] 
for collisions of a monochromatic 1 GeV neutrino beam 
on '^O. Instead of a sharp line at the incident energy, the 
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Figure 1: (color online) Reconstructed neutrino-energy distribution 
and its main components for a 1 GeV incident neutrino beam and '^0 
target. 

plot shows a broad peak due to Fermi motion. More im- 
portantly, there is a peak of non-CCQE events for which 
Ey reconstruction according to Eq. ([TJ is wrong. Many 
of them can be excluded by detecting the emitted pi- 
ons, but even with ideally perfect detection efficiency, a 
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remaining (CCQE-like) fraction of events correspond- 
ing to absorbed pions introduces a systematic error in 
the neutrino energy reconstruction that affects the de- 
termination of oscillation parameters (see for instance 
Ref. 121). Fake CCQE events can only be removed using 
a model dependent MC simulation. Therefore, a good 
insight into the dynamics of neutrino-nucleus (vA) col- 
lisions can be hardly underestimated. 

The importance of weak pion production (:7rP) for os- 
cillations studies is not limited to the contamination of 
CCQE samples. Neutral current (NC) jt" production 
(both incoherent and coherent) is a large background 
for Ve appearance searches. When one of the two pho- 
tons from a tt" decay escapes detection, the tt" cannot 
be distinguished from an electron born in a Vg induced 
CC interaction. Therefore, a precise determination of 
and the potential discovery of CP violation in the 
lepton sector requires that this background is reliably 
subtracted. Another Vg appearance background is pho- 
ton emission in NC, one of the relevant processes at low 
reconstructed neutrino energies where MiniBooNE has 
found an unexpected excess of events |[3l. 

The experimental information on neutrino interac- 
tions is rapidly growing. Detailed quasielastic (QE) 
scattering measurements have been published by Mini- 
BooNE at £v ~ 1 GeV (CCQE [T], NCQE 131) and NO- 
MAD at high (3-100 GeV) energies L6J (CCQE). Mini- 
BooNE has also reported single pion production cross 
section measurements on different channels: CCtt^ Q, 
CC7r° |81 and NCtt" |9|, including first single- and 
double-differential distributions. Special attention has 
been paid to coherent nP in both the CC and NC chan- 
nels by MiniBooNE [? !, SciBooNE GOlin] and NO- 
MAD lfT2l . These new data represent a major improve- 
ment and are challenging our understanding of vA inter- 
action physics. 

2. Modelling quasielastic scattering 

At the nucleon level, the weak interaction is defined 
by the current 



7aFl + 0-a(i(fF2 



qa 



-yajsFA ysFp 



UN , 



(2) 



written in terms of form factors (FF) Fi^2,a,p(Q^)- 
Vector-current conservation and isospin symmetry im- 
ply that Fi2 are given in terms of the nucleon elec- 
tromagnetic FF, extracted from electron scattering 
data |fT3l . Partial conservation of the axial current 



(PCAC) allows to relate Fp to Fa, usually parametrized 
in a dipole form 



Fa(Q^) = gA 



1 + 



-2 



(3) 



in analogy to the electric FF of the proton at low Q^. 
Once gA is fixed from neutron /3 decay, the axial mass 
Ma, related to the axial mean square radius (r^) - 
12/M^, remains the only unknown nucleon property in 
Eq. The value of Ma extracted from early CCQE 
measurements on deuterium targets is Ma - 1.016 ± 
0.026 GeV 113]. While one might be tempted to dis- 
trust this result based on experiments with low statistics 
and poorly know neutrino fluxes, there are good reasons 
to think that, at least at low Q^, Ma ~ 1 GeV. Indeed, 
there is a low energy theorem that relates n electropro- 
duction amplitudes to Fa at threshold and in the chiral 
limit. Using models to connect this theorem with data 
it has been found that Ma = 1.069 + 0.016 GeV [il4l . 
Moreover, applying a hadronic correction that can be 
precisely calculated at low using chiral perturbation 
theory |15 |, the resulting Ma - 1.014 + 0.016 GeV is 
even closer to the one from vd experiments. 

The vast experience acquired in electron-nucleus 
scattering studies has been applied to the vA case. The 
simplest model in the QE region, present in most event 
generators used in the analysis of neutrino experiments, 
is the relativistic global Fermi gas (REG) [16|. It as- 
sumes the impulse approximation (lA) according to 
which the interaction takes place on single nucleons 
whose contributions are summed incoherently. The 
struck nucleons have momentum distributions charac- 
terized by a Fermi momentum pf, and a constant bind- 
ing energy eg. Outgoing nucleons cannot go into oc- 
cupied states (Pauli blocking). Such a simple picture 
explains qualitatively inclusive QE electron scattering 
data but fails in the details. A better description re- 
quires a more realistic treatment of nuclear dynamics. 
Interacting nucleons do not have a well defined disper- 
sion relation but become broad states characterized by 
spectral functions (SF) 



S h,p — 



ImS 



7T (p2 -ml- ReE)2 + (ImE)2 ' 



(4) 



where Z is the nucleon selfenergy. The particle SF 
S p accounts for the interaction of the outgoing nucleon 
with the medium. The hole SF S h includes an 80- 
90% contribution from single-particle states while the 
rest of the nucleons participate in nucleon-nucleon in- 
teractions (correlations) and are located at high momen- 
tum IT71 [TSl ). In several models these correlations are 
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neglected. An example is the relativistic mean field 
(RMF) model lfT9l where the initial nucleons are treated 
as single-particle bound states whose wave functions 
are solutions of the Dirac equation with a cr-co mean 
field potential. The same energy independent real po- 
tential and Dirac equation are used to obtain the dis- 
torted wave function of the outgoing nucleons. The fact 
that the potential is real makes the model suitable for 
inclusive processes, in contrast to those with complex 
optical potentials more appropriate for exclusive reac- 
tions like A(v/, rp) because the imaginary (absorptive) 
part of the potential accounts for the flux loss towards 
other channels (see for instance Figs. 1-2 of Ref. EOl ). 

Another alternative to the RFG is the so called Lo- 
cal Fermi Gas (LFG) where the Fermi momentum de- 
pends on the coordinate through the nuclear density pro- 
file pf(r) = [(3 /2)7r^ p(r)]^'\ The LFG description in- 
troduces space-momentum correlations, absent in the 
GFG, that render more realistic the nucleon momen- 
tum distributions (see Fig. 6 of Ref. fTl]). A great 
advantage of LFG is that, owing to its simplicity, mi- 
croscopic many-body effects such as SF ||2T] |22l and 
long range random phase approximation (RPA) correla- 
tions Il23ll22li24 j are tractable in a realistic manner. 

2.1. The CCQE puzzle 

A common feature of all known calculations of the 
CCQE integrated cross section on '^C applying the dif- 
ferent theoretical techniques outlined above is that they 
underestimate recent MiniBooNE data. This surprising 
situation is illustrated in Fig. |2]for some of the model 
calculations collected in Ref. |25l and also the RFG. 
Theoretical results from quite different models lie on 
a rather narrow band (narrower than the experimental 
errorbars) clearly below the data: at - 0.8 GeV, 
cr,h -4.5-5 while cr^.,^ ~ 7 x 10"^^ cm^. The dif- 
ferences in the nucleon FF adopted in these calculations 
are minor In particular all take Ma ~ 1 GeV. 

Several interpretations of this discrepancy are cur- 
rently under debate. One points at the difficulty that 
the neutrino-flux determination represents and the pos- 
sibility that its absolute normalization has been under- 
estimated. On the other hand, according to the Mini- 
BooNE collaboration the systematic errors in the flux 
estimation have been determined by varying parame- 
ters within their uncertainties and accounting for cor- 
relations f29] so it is legitimate to expect that the er- 
rorbars account for the uncertainties in the flux nor- 
malization. Another strategy is to extract Ma from 
MiniBooNE data. In Ref. p|, a fit to the shape of 
the reconstructed distribution with the RFG model 
yielded Ma = 1.35 ± 0.17 GeV, which is much higher 



CCQE on C 




E, [GeV] 

Figure 2: (color online) Summary of CCQE total cross sections. Solid 
lines denote the models from Rets 1 18|, |26|, 127), (28||2TJ, (SI, (H 
and 1 22 1 in this order, as reported in Ref. |25|. The dash-dotted and 
dotted lines are RFG calculations with pf = 220 MeV, eg = 34 MeV 
and Ma = 1 and 1.35 GeV respectively. The dashed line is the result 
of Ref. 124J after adding the 2p - 2h contributions. The data points 
are from MiniBooNE (4). 



than the world average and the recent NOMAD result 
Ma = 1.05 ± 0.02( Stat) + 0.06( syst) GeV |6l at high 
energies. The integrated cross section computed with 
the new value of Ma is consistent with the normalized 
data as can be seen in Fig. |2] (dotted line). A similar 
fit but using a more elaborated distorted wave lA 
model, like the RMF sketched above but taking also 
into account nuclear correlations for the initial nucle- 
ons and with a different (real) potential for the outgoing 
ones, obtained Ma - 1.37 GeV |301. A better descrip- 
tion of the low region compared to the RFG was 
also achieved. With a state-of-the-art SF, the best 
fit and a good description of muon energy spectrum and 
angular distribution were found for an Ma as large as 
1.6 GeV |31 1. With a similar SF but fitting direcdy the 
measured flux-averaged double differential cross sec- 
tions {d^cr/dEfjd cos Ojj) taking into account the flux 
uncertainty and introducing a three -momentum cut of 
500 MeV to exclude the lA breakdown region, it was 
obtained that Ma = 1.343 + 0.060 GeV, lower but stifl 
incompatible with earlier determinations. 

A third possibility has been put forward by Martini et 
al. f24] . They have studied inclusive vA scattering in a 
LFG using RPA and taking into account two-particle- 
two-hole (2/7 - 2h) contributions, in particular some 
terms that are not part of the SF (see diagrams 2, 3, 3' 
in Fig. 1 of Ref. Jill). As shown in Fig.|2] with Ip-lh 
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excitations alone, the prediction of this model is consis- 
tent with the rest but the 2p-2h component turns out to 
be substantially large and capable of explaining the size 
of the cross section measured by MiniBooNE. An in- 
teresting prediction of the model concerns the antineu- 
trino CCQE reaction L32J : the different interaction pat- 
tern implies that, contrary to the neutrino case, 2p -2h 
excitations play a minor role. These suggestions should 
be further investigated by comparing the distribution 
and double-differential cross sections with data. A de- 
tailed validation with inclusive electron scattering data 
is also desirable. The role of meson-exchange currents 
(MEC) and relativistic effects needs to be elucidated. 
Work in this direction has already started with the cal- 
culation of vector MEC in the 2p - 2h sector with the 
RFGlOa. 



transition current can be cast as 



3. Modelling incoherent pion production 

Pion production in nuclei is incoherent when the fi- 
nal nucleus is excited (v;A — » InX). The first step to- 
wards a good description of nY' on nuclear targets is a re- 
alistic model of the elementary reaction (on nucleons). 
The most popular model for this process in neutrino in- 
teraction simulations was developed by Rein and Seh- 
gal |34|. It assumes that nP on the nucleon is dominated 
by baryon resonance excitation, which is described us- 
ing the relativistic quark model of Feynman, Kissinger 
and Ravndal [35] for resonances with invariant masses 
up to 2 GeV. The model originally neglected final lepton 
masses. This is a bad approximation at low but finite 
mass corrections in kinematics and currents have been 
recently investigated 1361 l37l [38l . More worrying is the 
poor description of electron scattering on the proton (see 
Fig. 2 of Ref. L39J and Fig. 2 of Ref [40 1) due to the use 
of unrealistic vector FF. On the other side, the wealth of 
pion photo- and electro-production data available from 
several experiments at MIT/Bates, MAMI/Mainz and 
specially JLab have been used to extract the electro- 
magnetic transition helicity amplitudes PTl l42ll . This 
valuable empirical information should be incorporated 
to the analysis of neutrino experiments. 

In contrast, there is almost no information about the 
axial part of the weak nucleon-to-resonance transition 
current. PCAC and pion-pole dominance of the pseu- 
doscalar FF can be applied to relate the axial coupling 
for the dominant contribution at low to the resonance 
nN decay coupling [off-diagonal Goldberger-Treiman 
(GT) relation, see for instance Appendix C of Ref. [21]]. 
In the few-GeV region, weak ttP is dominated by the 
excitation of the A(1232) resonance, for which the axial 



— ig"'i-q"y')+%g"'q-p' 
niN mi, 



^6 



■q"p"') +C^g»^' + ^q»cf 



75 uip) ■ 
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For small Q^, only the axial FF is relevant and some 
effort has been devoted to its extraction from ANL and 
BNL bubble chamber data. In Ref. ||43l C^(0) was ex- 
tracted from the ratio of the inelastic V/^d — » p^ pn 
and quasielastic Vf,d — » p^ p p distributions mea- 
sured at BNL. In the ratio neutrino flux uncertainties 
largely cancel out. The result C^iO) = 1.22 + 0.06 
is compatible with the GT value of 1.2. Assuming a 
simple dipole parametrization like in Eq. (3) for C^, 
Graczyk et al. El obtained C^(0) = 1.19 ± 0.08 and 
MA^ = 0.94 + 0.03 GeV by directly fitting daldQ^ 
for v^d p^ p n ANL and BNL data taking into 
account normalization uncertainties. Both studies in- 
cluded deuteron effects but neglected nonresonant back- 
grounds. The nonresonant contribution close to thresh- 
old is fully determined by chiral symmetry fA5\. Its in- 
clusion required a considerable reduction of Cg (0) - 
0.867 + 0.075 (with M^a = 0.985 ± 0.082 GeV) to 
describe the data, but the fit was done to ANL data 
alone, which are systematically below BNL ones. A re- 
cent reanalysis including also low-energy BNL data and 
deuteron effects finds a higher Cj (0) = 1 .00 ± 0. 1 1 , with 



Mak = 0.93 + 0.07 GeV [gl, but stifl implies a 20 % 
reduction of the GT relation. This result is in agreement 
with (0) = 0.96 obtained in the dynamical model for 
weak ttP of Ref. [47 [ after a 30 % renormalization of the 
constituent quark model prediction by meson clouds. 

When ttP takes place inside the nucleus, the strong- 
interacting environment leaves a big imprint on the ob- 
servables. First of all, the elementary amplitude is mod- 
ified in the medium by the presence of the nuclear mean 
field and, most importantly, due to the modification of 
the A(1232) resonance, whose mass gets shifted and 
width increased because of absorption processes involv- 
ing one (AA^ — > N N) or more nucleons. In addition, 
the produced pion can be absorbed or scatter with the 
nucleons with and without charge exchange. In the 
few-GeV region a large number of states can be ex- 
cited so that the only feasible way of describing the ex- 
clusive final system is with a semiclassical treatment. 
The most common framework to deal with this is an 
intranuclear cascade but transport theory has also been 
used. Reviews and short descriptions of the different 
MC models applied to vA interactions can be found in 
Refs. 1481 |49l . Here I shall only refer the main fea- 
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tures of nuclear nP calculations undertaken by theory 
groups. In the NuWRO event generator |50| pions are 
produced via A(1232) excitation. Heavier resonances 
are only accounted for in the duality-inspired nonres- 
onant background. Pion propagation is accomplished 
with an intranuclear cascade, with scattering probabil- 
ities determined by nN vacuum cross sections. Pion 
absorption is fixed according to pion nuclear absorp- 
tion data. The model of Ahmad et al. 151] takes also 
into account the A(1232) in-medium change only in the 
production amplitude while the pion cascade uses vac- 
uum cross sections. Finally, the Giessen Boltzmann- 
Uehling-Uhlenberg (GiBUU) model is a semiclassical 
transport model in coupled channels successfully ap- 
plied to photo-, electro- and hadron-nucleus reactions 
and recently extended to vA collisions Il28l . Not only 
the A(1232) but all the baryon resonances with masses 
up to 2 GeV can be weakly excited and are explicitly 
propagated. The medium modifications equally affect 
the production mechanism and secondary interactions. 

3.1. The incoherent pion production puzzle 

MiniBooNE has measured the ratio 
cr(CC7r+)/cr(CCQE-like) on CH2 as a function 
of Ey |,52J. There is an uncertainty in the neutrino- 
energy reconstruction but, on the other side, this 
observable does not depend on the neutrino flux 
normalization. The three models described above have 
been employed to calculate this ratio: see Fig. 8 of 
Ref. f44|. Fig. 1 of Ref. ||53l and Fig. 2 of Ref. 1541 . 
The comparison shows a good agreement at the lowest 
energies that gets progressively worse as the energy 
increases, with NuWRO |44| getting better agree- 
ment while the GiBUU calculation |54| exhibits the 
largest discrepancies. Recall however that the three 
calculations have used M/^ ~ 1 GeV for the CCQE 
cross section that enters the denominator, which is 
clearly insufficient to explain MiniBooNE CCQE data 
as shown in Fig. [2] Therefore it turns out all these 
models underestimate considerably MiniBooNE ;rP 
cross section. Indeed, the comparison of the NCtt" pion 
momentum distribution obtained by GiBUU f55l with 
the corresponding MiniBooNE data |9J reveals a good 
agreement in the shape but a factor of ~ 2 discrepancy 
in the normalization which is a major challenge in our 
understanding of weak ttP. 

4. Modelling coherent pion production 

Pion production in nuclei is coherent when the nu- 
cleus remains in its ground state (v/A — > In A, with 
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Figure 3: Pion momentum distribution for NC;r production on CH2 
averaged over the MiniBooNE flux |29 | calculated with the GiBUU 
transport model 1551 . Data are from Ref. 0. 



JT - in the CC case and n - tt" in the NC one). It 
takes place predominantly when the momentum trans- 
fered to the the nucleus is small so that the latter is less 
likely to break. Weak coherent ttP has a very small cross 
section compared to the incoherent process, but rela- 
tively larger than equivalent reactions induced by pho- 
tons or electrons due to the non vanishing contribution 
of the axial current at the relevant kinematics ll56l . 

The pioneering model of Rein and Sehgal (RS) ISTl 
uses PCAC in the = limit to relate neutrino- 
induced coherent ;rP to pion-nucleus elastic scattering, 
given in terms of the pion-nucleon cross section. By 
taking the - limit, the RS model neglects im- 
portant angular dependence at low energies lf58l . This, 
together with the fact that the description of the pion- 
nucleus elastic cross section is not realistic (see Fig. 2 
of Ref. 1581 ). results in cross sections well above the 
experimental data. An alternative approach based on 
PCAC directly uses the experimental pion-nucleus elas- 
tic cross section ||59l l60l . In this way the treatment 
of the outgoing pion is improved, but a spurious ini- 
tial pion distortion, present in pion-nucleus elastic scat- 
tering but not in coherent ttP, is introduced. With this 
method the cross sections are smaller and more com- 
patible with the experiments. As the energy increases, 
the deficiencies become less relevant and the RS, and 
PCAC models in general, becomes realistic. In fact, 
the recent NOMAD measurement of the NCtt** coher- 
ent cross section at high (3-100 GeV) energies found 
cr = [72.6 + S.l(stat) ± 6.9(syst)] x lO--*" cm^ which is 
consistent with the RS prediction of 78 x 1 0"*° cm^ |Il2l . 

Microscopic approaches meant to work in the A re- 
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gion have also been developed. They combine the A 
excitation picture of weak ;rP on the nucleon or the 
more complete models of Refs. Il45l l47ll with the A- 
hole model of pion-nucleus interaction. As the nu- 
cleus remains in its ground state, a quantum treat- 
ment of pion distortion is feasible by means of the 
eikonal approximation 11611 . the Klein-Gordon Il62ll56l 
or the Lippmann-Schwinger f6T| equations with real- 
istic pion-nucleus optical potentials. In analogy to the 
incoherent case, it is found that medium effects and 
pion distortion reduce considerably the cross section 
and shift the peak to lower pion-momenta (see for ex- 
ample Fig. 2 of Ref. Il56l ). Nonlocalities in the A 
propagation were found relevant |64| and have been 
explicitly taken into account in Ref. 1.6 3 J but not in 
Refs. ||6ll|62l|65l|56l|66l, where they are (at least par- 
tially) included via the empirical A mass shift f63','66l. 

The size of the coherent ttP cross section is very sen- 
sitive to the value of (0) 1651 . This can be easily un- 
derstood from the fact that at = (not reachable in 
the CC case with nonzero lepton masses), where most 
of the strength of this reaction is concentrated the only 
FF that contributes is . Therefore at low energies, 
to a good approximation, the cross section scales with 
[Cj (0)]^. For this reason, establishing the degree of vi- 
olation of the off-diagonal GT relation for the N-A tran- 
sition as discussed in the previous section is crucial to 
predict the coherent pion production strength at low en- 
ergies. 

4.1. The coherent pion production puzzle 

SciBooNE has recently published a measurement of 
the CC/NC coherent pion production ratio employing a 
neutrino beam of {e}) ~ 0.8 GeV HI]- They find that 



cr(CCcoh7r+) ^ ^^_,o 30 
cr(NCcoh;rO) ' ""'^^ ' 



(6) 



The theoretical models outlined above predict different 
cross sections. Different are also the predictions from 
various MC generators even if all implement the RS 
model; see Fig. 8 of Ref. ||25]| for a compilation. How- 
ever, the ratio is always in the CC/NC ~ 1 - 2 inter- 
val for Ey ~ 0.8 GeV, which is in strong contradiction 
with the SciBooNE result. Indeed, neglecting final lep- 
ton masses and in the = limit this ratio is equal 
to 2; in a realistic situation one expects only some cor- 
rections to this result. For example, in Ref ||66l the 
ratio of the flux-averaged cross sections is found to be 
1.46 + 0.03, below 2 as one would expect, but far from 
the SciBooNE value. 



5. Modelling single photon emission 

Like weak nP, NC single photon production can be 
incoherent (vA — > vyX) or coherent (vA — » vyA). 
These processes were studied in Ref. Ii67i with effec- 
tive Lagrangians including nonresonant Compton-like 
diagrams (Fig. 1 of Ref. |67|), f-channel contributions 
from the anomalous vertices with tt, p and a> meson 
exchange (Fig. 2 of Ref. ll67l ) and A(1232) excitation 
terms (Fig. 3 of Ref. [67l). The cross section on sin- 
gle nucleons and the coherent nuclear one were com- 
puted. The A(1232) contribution was found to be the 
largest one in both coherent and incoherent cases. The 
f-channel part is dominated by oj exchange. It was 
found that MiniBooNE's strategy of normalizing A — > 
Ny events to the measured A ^ Nn rate is a reason- 
able first approximation for the incoherent reaction but 
not for the coherent one. Within the uncertainties, the 
model seems capable of providing enough photons to 
cover the low-energy excess found by MiniBooNE 1681 
but a precise account of nuclear effects is necessary. 

6. Concluding remarks 

New neutrino interaction data with high statistics, ac- 
companied by a better understanding of the neutrino 
fluxes are becoming available from several experiments. 
The comparison with theory for some of the reaction 
channels that are relevant for oscillations studies and 
interesting from the perspective of hadronic and nu- 
clear physics reveals discrepancies that await explana- 
tion: fitting the new results with the available param- 
eters is a dangerous strategy in the long term. The 
general tendency is that theory underestimates data but 
some flux-normalization independent quantities like the 
cr(CC7r^)/cr(CCQE - like) ratio are also not well de- 
scribed by state-of-the-art models. CC coherent tt^ pro- 
duction turned out to be elusive even if all theoretical 
models predict its cross section to be larger than the NC 
coherent tt" one. Future data and ongoing theoretical 
work shall be quite helpful to clarify the situation. In 
order to achieve the precision goals in neutrino oscil- 
lation measurements and to reliably extract information 
about the axial properties of the nucleon and baryon res- 
onances it is crucial that the current theoretical develop- 
ments are implemented in the event generators used in 
the experimental data analysis. 
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